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Abstract 

Recent results on the stationary state Fluctuation Theorems for work and heat fluctuations 
of Langevin systems are presented. The relevance of finite time corrections in understanding 
experimental and simulation results is explained in the context of an exactly solvable model, 
namely a Brownian particle in a harmonic potential, which is dragged through the surrounding 
fluid. In this model, work fluctuations obey the conventional form of the fluctuation theorem while 
heat fluctuations satisfy an extended form. The connection with other work in recent literature 
is pointed out, and further generalizations are suggested. 



1 Introduction 

In the last twelve years, the field of non-equilibrium statistical mechanics has been given a boost 
by the discovery of a number of new, general results regarding the fluctuations of thermodynamic 
quantities such as heat, work and entropy. This paper is not intended as a survey article, rather, we 
will restrict ourselves to non-equilibrium stationary state Fluctuation Theorems (FTs). A stationary 
state FT was first discovered in [1]. This FT is a quantitative relation between the probability to 
observe a positive entropy production a versus a negative one —a over a time interval r. Given the 
probability distribution 7r r of a fluctuating quantity p, scaled such that its average is equal to one, the 
Conventional Fluctuation Theorem (CFT) (also called the Gallavotti-Cohen Fluctuation Theorem) 
takes the form [2] 

1 . T T (P) , * /-.n 

hm — log — = a + p tor p < p , (lj 

r^oo r 7T T (— p) 

where a + is the average (dimensionless) entropy production rate over positive times, p = a/a + , and 
p* is an upper bound on the values of p for which the CFT holds. 

Whether the CFT holds depends on the kind of system and on what the quantity p represents. 
Choices for p that have been used in the literature include the phase space contraction rate [2], a 
"dissipation function" [3], an "action functional" [4], and, more concretely, work, heat and entropy 
production. The value of p* can also depend on the quantity under consideration. The stationary 
state Fluctuation Theorem holds quite generally (i.e. for a large class of deterministic [1-3] and 
stochastic systems [4,5]), and has been verified in laboratory experiments [6,7]. 

The behavior for values p > p* has recently been elucidated for specific systems, and led to an 
Extended Fluctuation Theorem (EFT) for heat [8-10], which we will discuss in detail below. 

In this paper, we will give an overview of some interesting results of the stationary state Fluctu- 
ation Theorem for an exactly solvable model of a Brownian particle in a harmonic potential which 
is dragged through a surrounding fluid. While this model is based on a Langevin equation, the 



conclusions from this model also pertain to certain purely deterministic systems, as we will point 
out. 

We will first introduce the model in Sec. [2] and discuss its exact solution in Sec. El We will 
then turn our attention to thermodynamic quantities of the dissipated heat and the supplied work, 
and their fluctuations, both for finite and infinite time in Sec. HI In Sec. [5] we will show that this 
results in two kinds of stationary state fluctuation relations: the Conventional Fluctuation Theorem 
(CFT), for work fluctuations, and an Extended Fluctuation Theorem (EFT), for heat fluctuations. 
The paper ends with a discussion in Sec. [6] of the results and generalizations that have appeared in 
the literature. 



2 The model 



Consider a spherical Brownian particle of mass m and diameter D suspended in an equilibrium fluid 
with viscosity rj and temperature T. The motion of this Brownian particle can be described by the 
Langevin equation 



mx = —av + £, 



(2) 



where x is the position of the particle, v = dx/dt its velocity and a = 3irr]D the Stokes friction 
coefficient. Furthermore, £ = is a Gaussian fluctuating force whose average (over realizations) 
(£(t)) = and whose time correlations are delta-correlated, i.e. 



<£(*i)£(*2)> =2k B TaS(t 1 -t 2 )l, 



(3) 



where &b is Boltzmann's constant and 1 is the identity matrix. 

Consider now the situation in which an external force acts on the Brownian particle (but not on 
the fluid). Under the assumption that the fluid is unperturbed such that a and the properties of £ 
are not changed, the Langevin equation gets modified to 



mx 



F(x, t) — av + 



(4) 



where F(x,t) is the external force, which may depend on position and time. 

Of particular interest here is the case in which the force derives from a moving harmonic potential 
U, i.e., 



F(x,t) 



dU(x,t) 
dx 



n(x — x*(t)), 



(5) 



where the potential U(x,t) — Ki | X X *(i)| 2 /2 and the center x*(t) of the potential changes in time 
according to 



x*(t) = v*t. 



(6) 



This corresponds to the particle being trapped in a harmonic well, the center x* of which moves 
with a constant velocity v* , to which we will restrict ourselves here. This closely resembles the 
experimental system of a Brownian particle in a moving optical trap studied by Wang et al. [11]. 

One also needs to specify the initial condition of the system. In time, the distribution of x and v 
approach a non-equilibrium stationary state distribution _f ness asymptotically. It is not hard to show 
that 



JnessKX, V) = — exp 

Z7T 



v*\ 2 + k\x 



v*(t - 



(7) 
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where (3 = 1/(A;bT). Note that if v* = this reduces to the equilibrium distribution. To focus 
on the non-equilibrium stationary state fluctuations, one can take as an initial condition that the 
distribution at time t = 0, denoted by /o, is equal to / neS s- 

This system exhibits many features of general non-equilibrium systems, namely the presence of 
a current, energy dissipation, and entropy and heat production due to the continuous input of work 
from the outside. Thus it can serve as a prototypical non-equilibrium system. The advantage of 
looking at this particular system is that it can be solved exactly (as was first exploited for work 
fluctuations by Mazonka and Jarzynski [12]). In fact, this system can also be solved exactly for a 
general the time dependence of the center x*(t) of the harmonic potential, as was shown in [13]. 

3 The exact solution 

The exact solution of the Brownian particle system introduced above proceeds through a time- 
dependent change of coordinates that brings the equations into a form in which the force does not 
depend on time, and whose solution is known. 

As explained in [13], to eliminate the time dependence of the force, such a coordinate transfor- 
mation should take the form 



where x is any solution of the equations of motion in the absence of the fluctuating force. One such 
solution is 



Note that the new coordinate system corresponds to a frame that is co-moving with the fictitious 
noiseless trajectory x(t), and that along this fictitious trajectory, the friction force —av* is exactly 
balanced by the harmonic force —k(x — x*) = av* . 

In the co-moving frame, the Langevin equation becomes [cf. Eqs. — (JSJ)] 



where V = dX jdt = v — v*. Equation (|10p is of the same form as Eq. (jj]) but with v* = 0, so it 
describes the system with a stationary harmonic potential. This system is well-studied and is an 
example of a so-called Ornstein-Uhlenbeck process [14]. Roughly speaking, an Ornstein-Uhlenbeck 
process is a stationary Markov process that can be described by a Langevin equation with linear 
forces and a Gaussian random force term. Equation (jlOH is Markovian because the time-derivatives 
of X only depend on quantities at that same time. As a consequence of the Markovian nature of the 
process, all its properties can be expressed in terms of the Green's function Gt(X' , V; X, V) which 
is the probability, in the co-moving frame, for the particle to be at position X' with velocity V given 
that it was at position X with velocity V a time t earlier. For Ornstein-Uhlenbeck processes, the 
Green's function can be solved exactly and is Gaussian in both the initial coordinates (X, V) and 
the final coordinates (X' ,V). It is this Gaussian property that enables one to obtain many explicit 
results, as we will see. 

Although this system can be solved exactly, to avoid unnecessarily complicated expressions, we 
will consider here only the overdamped case, which means taking the limit m <C a 2 /k. In Eq. (|10[) . 
this limit can be taken by setting m equal to zero, so that it reduces to = —kX — aV + or 



X(t) = x(t) - x(t), 



(8) 




(9) 



mX = —kX — aV + £ 



(10) 



X = -—X + a' l $, 
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where r r = a/re is the (average) relaxation time of this Langevin equation. Like the general case, 
this equation also describes an Ornstein-Uhlenbeck process, whose Green's function is given by [14] 



/?re 



,1/2 



/3re \X' 



G t (X';X) 



2tt(1 - e" 2 */^) 



exp 



2 



I _ e -2t/r r 



(12) 



Given a set of times < t\ < < ■ ■ ■ < t m , one can write for the joint probability of the particle 
to be at position X n at time t n (with to = 0): 



Because each Green's function is Gaussian in both the initial and final positions [cf. Eq. (I12p ]. and 
the initial distribution in Eq. ([7]) is also Gaussian in the initial position Xq, the distribution P is 
also Gaussian. 

4 Thermodynamic fluctuations 

4.1 Definitions 

While on a macroscopic level, thermodynamic properties such as work, heat and entropy are well- 
defined, this is not the case on a mesoscopic level. To understand that there is an ambiguity, consider 
trying to determine the "energy of the system" for the Brownian particle under consideration. Should 
this include the potential energy of the external harmonic force, or should this only be the (kinetic) 
energy of the Brownian particle itself? Or consider the "work done on the system" . Should this be 
the work done by the harmonic force on the particle, or should this be the work needed to move the 
harmonic potential? The first choice is closely related to the mechanical definition of work, but it 
leads to work being done even in equilibrium (although its average would be zero), a problem that is 
avoided by using the second choice. There are even more ambiguities when studying dissipation in 
non-equilibrium systems. When should energy be regarded as dissipated? It would be ambiguous to 
say that the dissipated energy is that part of the energy input that cannot be used to perform work. 
For instance, the harmonic potential energy can be used to store energy by displacing the Brownian 
particle from its origin. Is this dissipated energy? A clever experimentalist may be able to extract 
useful work from this energy, but should the definition of dissipated energy depend on the ability of 
an experimentalist? 

The ambiguity is thus due to what to call the system, and what forms of energy to consider as 
inaccessible or dissipated. To avoid this ambiguity, we will now define what we will mean by "work" 
and "heat" in the rest of the paper. The work will be defined as the work that is required to keep the 
harmonic potential moving. Because at time t the harmonic potential exerts a force —n(x{t) — x*(t)) 
on the particle, the particle will exert a reaction force n(x(t) — x*(t)) on the harmonic potential 
(or on whatever device generated it). To keep the potential moving therefore requires (at least) a 
compensating force — n(x(t) — x* (t)) to be applied to the potential. Since the velocity of the potential 
is x* = v* (cf. Eq jSD), the mechanical work done by this compensating force during a time interval 
r is 



Note that this definition of "work on the system" uses the particle plus the harmonic potential as 
the system. 

On the other hand, we will define the heat as the energy which is given off to the fluid, taking the 
point of view that since the energy of the fluid is not explicitly present in the Langevin description, 



m 



P(X , X U ...,X m )= fo(X ) J] Gu-t^ {Xi- X;_!). 



(13) 



8=1 




(14) 
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we can suppose it to be unaccessible or dissipated energy. This leads to the following expression for 
the heat: 



Q T = W T - AU T (15) 

where AU T = U(x(t)) — U(x(0)) is the change in harmonic potential from time to time r. In 
words, this states that work can be converted both into potential energy of the Brownian particle, 
or into heat dissipated into the fluid. Note that in the general (non-overdamped) case, the change in 
kinetic energy of the particle should also be subtracted from W T to obtain the heat. 

4.2 Fourier transforms 

Given the definitions of the work and the heat and the Green's function of the process, we can study 
the probability distributions of their fluctuations. Once the probability distributions are obtained, 
it can be checked whether they obey a stationary state Fluctuation Theorem. In addition (and in 
contrast to the majority of theoretical works on the subject), the explicit solution available for this 
system allows us to consider finite time corrections to these theorems as well. These corrections 
give information about the time scales on which these stationary state Fluctuation Theorems can be 
observed. 

An essential ingredient in the calculation will be that any linear function of a Gaussian distributed 
quantity is again a Gaussian distributed quantity, so that the distribution of such a quantity is fully 
determined by its mean and its variance. While the work W T in Eq. (|14p is a linear function of x(t), 
making it Gaussian distributed, the heat Q T in Eq. (|15p is not a linear function of x(t), but is in 
fact a quadratic function of x(t) and consequently not Gaussian distributed. This quadratic nature 
prevents a straightforward calculation of the heat probability distribution function, but still allows 
the computation of its Fourier transform, using the following trick. 

Consider first the joint probability P*(W T , £c(0), x(t)) which gives the probability that the particle 
moves from position x(0) to position x(t) and an amount W T of work is done in the process. Since 
these three quantities are all linear functions of x(t), their joint distribution is Gaussian, i.e. of the 
form 

-^(y-MyV-^iy-M) 

™ = <16) 

where for brevity we denoted y = (W t ,x(0),x(t)), M as the mean of y and V as the variance 
matrix of y. 

Using P* , the distributions of W T and Q T can be expressed as 

Pf(PT) = jdyP* T {y)5(W -W T ) (17) 
P?{Q) = JdyP;(y)5{Q-W T + K \x(T)-x*(T)\ 2 /2- K \x(0)-x*(0)\ 2 /2), (18) 
where x*(t) = v*t as given by Eq. ©. Taking the Fourier transform of these equations yields 

pW(q) = fdyP* T {y)e^ w - w ^ (19) 

pQ( g ) = /■ dyJ ^ (y)e iff(O-W- T +«|«»(r)-«*r|V2-*|«»(0)|V2) j ( 20 ) 

where q is the Fourier variable. The integrals on the right hand sides are all just Gaussian integrals, 
which can be performed explicitly once M and V are determined. The results for the stationary 
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state are then [9, 13]: 



P?<a) 



a+rq 



cr+Tq[i-q 



i— o 1 ; — 

1- 



r/rt-x 



r/T r (l + (l-e 



-2t/t,-^2) 



[1 + (l- e - 2 ^r)^2j3/2 



21 



(22) 



where g = ksTq is the dimensionless Fourier variable and <r+ = a/3|v*| 2 is the average dimensionless 
heat production rate. 

4.3 Saddle point approximations and singularities 

To obtain the actual probability distributions of work and heat requires the inversion of the Fourier 
transforms P™ in Eq. (HQ and P? in Eq. ([22]). Since Py is Gaussian, its inverse is easily obtained: 



(f)W-a + TY 



r/r r ' 



— 



2na+T 1 



l-e~ T / T r 
r/r r 



(23) 
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In fact, it is not necessary to use the Fourier transform to obtain this result for P^ , but it will be 
useful to have the Fourier transforms of both the work and the heat distribution to understand the 
origin of the peculiarities of the behavior of the heat fluctuations. 

The quantity of interest for the stationary state Fluctuation Theorems to be derived in Sec. 5 is 
not the probability of the work fluctuations but that of the scaled work fluctuations p = (3W/{a + r) 
(which compares the dimensionless W T , (3W T with the average dimensionless entropy production 
over a time r). Its probability distribution ir^ ' (p) is related to that of W by a Jacobian: tt^(p) = 
pW(W)dW/dp = pW(pw/(a + T))(a + T/P), so that 



T + T(p-ir 



tJtt 



*f(p) 



2vr 1 



\-e- T l T r 
t/t t 



(24) 



/(*+r) 



The Fourier inversion of P-? is more involved because no closed form is known for the Fourier 
inverse of Eq. (|22p . We are therefore left with the following formal expression for ir® {p) '■ 



7T«(p) 



£±I 

2vr 



a + Tq[i-q\ 



dq e 



-iqa+Tp 



r / r r 



r(l+(l-e- 2T / T '-) ( j2) 



[l + (l-e- 2 ^')g 2 ] 3 / 2 



(25) 



We are especially interested in the large r behaviour of this function. The presence then of a large 
parameter in the exponent in the integrand allows us to obtain an approximate result using the 
saddle point approximation. To exhibit the large parameter r in Eq. (|25p explicitly, we write 



a i r 



2vr 



7T«(p) 

where 



(26) 



2q\\ 



-t/t t 



(r/r r )(l + (l- e - 2 ^)5 2 ) 



2r 



log 



1 + (1 



-2t/t t \ -2 
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(27) 

In the saddle point method, the main contribution to the integral is found by considering the 
maximum q* of the function hr(q)- A complicating factor is that hr is complex valued. If one takes 
q* to be the point on the real line where hr has a maximum real part, then the integrand can be 
highly oscillatory there, causing the contribution of this otherwise dominant point to be washed out. 
A maximum that does not suffer from this oscillatory behavior can be found, but usually not on the 
real axis; indeed the solution of 

dh?(q*) , N 

„ W ' = 28 
dq* K ' 

lies here in the complex plane. 

An interesting property of complex functions is that they have no true maximum. According to 
the Cauchy-Riemann identities for complex functions, the point q* defined by (128j) has the property 
that there is a path through it called the path of steepest descent, on which has a constant 
imaginary part and its real part has a maximum, while along another (orthogonal) path the function 
has a minimum (also with constant imaginary part); thus q* is a saddle point of ■ The saddle 
point character of q* does not interfere with it being the dominant contribution to the integral, as 
long as we can deform the original line of integration (the real axis) into a contour which overlaps 
with the path of steepest descent through q* . Furthermore, it can be shown that along this line the 
function does not oscillate, so that the point q* truly gives the dominant contribution. Note that 
during the deformation of the integration contour, the endpoints (±oo on the real axis) must be 
fixed and no singularities of the function must be crossed. The value of integral along this deformed 
contour is then exactly equal to that of the original integral. To get the saddle point approximation, 
one expands the function hr{q) around q* up to second order in 5q = q — q*, and integrates along 
the deformed contour. As t increases, the function becomes more peaked around the saddle point, 
and the saddle point expression becomes an increasingly good approximation to the integral. 

To illustrate this method, let us first apply it to the work distribution, which is simpler than the 
heat distribution and for which the result can be compared with the exact result in Eq. (|24p . From 
Eq. (|2ip one then obtains 

«V(p) = ^f r dqJ*?®, (29) 



where 

/ 1 _ p-r/T r 



h^(q) = a+ { -iqp + q 



9 1 



r/r r 



(30) 



Since this is a quadratic function, the saddle point q* can be obtained analytically from dh^ /dq* = 0, 
with the result 

" ' ' ~ '' (31) 



2(1- l-e- T/Tr 



Note that the saddle point lies on the imaginary axis. Noting also that h!^(q) has no singularities as 
a function of q, one can deform the integration contour (i.e. the real axis) to go through q* . Around 
q* the function fi^ has the following behavior 

+ a « = — < iri (1 :5^ + a? ( i - I > ■ (32) 



1 For more general forms of h, there would be 0(Sq 3 ) correction terms. 
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Figure 1: Saddle points and paths of steepest de- 
scent of the functions (left) and (right) 
for a + = 1, r = 4r r and p = —2 (see text). Solid 
circles represent saddle points and solid curves 
S\, . . . ,Sjy drawn through these are the paths of 
steepest descent. For , shown on the left, only 
one saddle point B exists and the (dashed) real 
axis R can be deformed to the indicated path of 
steepest descent 5b going horizontally through 
B, without crossing a singularity. For hr shown 
on the right, wiggly lines represent branch cuts, 
two singularities appear at ±i, and there are four 
saddle points A, B, C and D. The path of steep- 
est descent Sb through B is the (only) way to 
deform the real axis to a path of steepest de- 
scent without crossing a singularity. 



Note that for real 5q, this expression is real (i.e., non-oscillatory) and has a maximum at 5q = 0: the 
path of steepest descent therefore crosses q* parallel to the real axis. Thus we can indeed deform the 
original line of integration (the real axis) by shifting the integration contour up or down (depending 
on p) until it crosses q*; this changes the end points ±oo, but the integrand vanishes there, so that 
the integral over the shifted line is still the same as for the original integration line. This is sketched 
in left panel of figure 1. The final integral now becomes 

_ <7-|_T(l- 

■p. v 

2tt 



^(P) = T-^-W d5 q e -° + \ —r. (33) 



After performing the Gaussian integral, one recovers Eq. ([24 

We now turn to the distribution of heat fluctuations. To find the saddle point in this case, the 
derivative of hr(q) has to be equal to zero. This derivative is rather complicated, in contrast to the 
simple linear function that it was for the work functions. To make the expressions in the following 
a bit more manageable, let us assume that r is so large that e~ T ^ Tr may be neglected in Eq. (I28p 
(by comparison with numerical methods, we found that this typically leads to accurate results for 
r > 3r r [9]). In that case, becomes 

h?(q) = a + {-iqp + q{i - q) + ^ - f) } - ± log(l + ?). (34) 

while 

= a + { -* + i - 2i + (t/t *£ - & + (t/t w } - -jjM-. (35) 

Note that the naive limits r — > 00 of h® in Eq. (1341) and in Eq. (|30p coincide. However, the terms 
in Eq. (|34p that are then neglected have singularities, which have to be taken into account, because 
even though the factor 1/r may tend to zero as r — » 00, it multiplies an expression that can, because 
of the singularities, go to infinity, making the limit r — > 00 ill-defined. 

The saddle point equation dh? /dq = leads, using Eq. (f35l) . to a fourth order polynomial equa- 
tion. Such an equation has four solutions; hence there are four saddle points, in contrast to the Fourier 
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transform of the work distribution, which had only one. Which saddle point (s) to use depends on 
how we can deform the real axis integration contour to include saddle points without crossing any 
singularities. As Eq. (|34|) shows, there are singularities at q = ±i. The point q = i is a branch point 
associated with the logarithm, while q = —i is in addition a simple pole. The structure of the saddle 
points is shown in the right panel of figure [TJ As a function of p, there is only one saddle point that 
we can reach by deforming the real axis without crossing singularities. This saddle point lies on the 
imaginary axis and is in a sense trapped between the singularities at ±z; this is quite different from 
the saddle point for the work fluctuations [i.e. p)/2 for large r, cf. Eq. (|3ip ] which is unbounded 
as a function of p. 

The full calculation of the Fourier inverse of iv would be too detailed for this short paper and 
can be found in [9]. In an asymptotic expansion for large r, the result is 



7T%) 



^b±M e -a +[TV -T P ] + 3/2 ifp<-l. 

OD7T r 



Wa+T e -"+(l-P) z [T+2(l-p)/(3-p)r r ]/4 



[(3-p)(l+p)F 



2 



if - 1< p < 3. (36) 



°+fjTr -°+ [(p-2)r-V8(p-3)rr r + 2^|2r r ] . f > ^ 



32tt 



Note that this function has a Gaussian center for —1 < p < 3 and exponential tails for p < — 1 and 
p > 3, i.e., for large r, one has ir T ~ expfo+rp], 7r T ~ exp[— a + r(p — 1) 2 /4] and 7r T ~ exp[— <r + r(p — 2)] 
for p < — 1, — 1 < p < 3 and p > 3, respectively. Qualitatively, this result can be understood as 
follows: if p is such that the saddle point q* is relatively far away from the singularities, the function 
and its relevant saddle point resemble those of the work fluctuations, i.e., q* scales linearly with p 
and the Fourier inverse is consequently Gaussian in p. Since for the work fluctuations q* = i(l—p)/2, 
the values of p for which q* is not close to the singularities ±i lie within the interval p £ [—1,3]. For 
values p > 3, the work saddle point lies below —i, but the heat saddle point has to stay within [—i, i]. 
As a result it is basically stuck just above —i, i.e., q* « —i, yielding, cf. Eq. (f25|) . an exponential 
dependence of ir?(p) on p. Likewise, for p < —1, the heat saddle point is stuck just below i and 
again an exponential tail of the distribution is present. 

Eq. (|36p also includes finite r effects, obtained using the saddle point approximation. Note that the 
approach to the large r behavior has a different character for different values of p. This corresponds 
to the different character of the singularities at dbi (i.e., pole vs. branch point). 



5 Stationary State Fluctuation Theorems 

Given the distribution functions, one can now investigate whether the stationary state Fluctuation 
Theorem holds. The quantity considered in the Fluctuation Theorems is the fluctuation function 
fr(p) = ^7 l°g[ 7r T(p)/ 7 i"r( — p)]- This function measures the asymmetry between positive and negative 

ir^(p) and Eq. ([21 



fluctuations in the quantity p. When p represents work, vr T (p) 
P 



gives: 



f7ip) 



i 



t/t t 



(37) 



As r — > oo, this tends to p from above (cf. figure 2), hence the Conventional Fluctuation Theorem 
(CFT)is satisfied by the work fluctuations. For these fluctuations, the characteristic time is r r . 

On the other hand, considering the heat fluctuations and using Eq. (|36p . one finds an Extended 
Fluctuation Theorem (EFT): 



f?(p) 



p- { ^ + 



3 In 



3-p*+2p 
3-p 2 -2p 



2 + 



8(2-3) 



r/r r 



T r 2(7+ T 

(p+l)(-2p 2 +8p-10) 
4(3-p)r/r r 

ln(r<T 2 Tr) 



rT t\ — 1 In (3-p) 3 (l+p) 3 (p-l) 
°+' ) 2 111 576 



2ct + t 



if < p < 1 
} if 1< p < 3 
if p > 3 
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Figure 2: Conventional (CFT) and Extended 
Fluctuation Theorem (EFT). W T= fi n it e and 
W T=ao represent the CFT and <3 T =fimte and 
Q T=OQ the EFT for finite and infinite times, re- 
spectively. The arrows indicate how the infinite 
time limit is approached, i.e., for the work fluc- 
tuations from above, and for the heat fluctua- 
tions from below for small p and from above for 
large p. 



(38) 

Note that the behavior for negative p can simply be found from the antisymmetry of the function 
f?(p), i.e. f?{-p) = -f?{p). 

To illustrate the difference between the work and heat fluctuation functions f}^ and fr , we 
sketched these functions as a function of p in figure [2 One sees that in the case r — ► oo, for p < 1, 
the fluctuation functions tend to a straight line with slope one. This behavior persists beyond p > 1 
for f^, while fr bends down away from the straight line and settles to a plateau value of 2 for p > 3. 

Thus, even in the limit r — > oo, Eq. (|38|) for heat fluctuations differs from the CFT for work 
fluctuations when p > 1, since the fluctuation function has for heat a plateau at a value of 2 rather 
than that it increases without bound as for work. Alternatively, one could say that fr does satisfy 
Eq. (pQ), but with p* = 1, whereas f\^ has p* = oo. We note that for the heat fluctuations there is in 
addition to r r a relaxation time 0+ 1 , which can exceed r r . 

The plateau value 2 in the EFT for p > 3 (cf. Eq. ([38]) ) is a consequence of the exponential tails 
of the distribution . These tails are related to the presence of singularities at ±z, which limit the 
position of the saddle point q* , such that for p > 3 one has q* « —i, while for p < —1, q* pa +i. As 
the saddle point gives the dominant contribution to the Fourier inverse, thus leads to the exponential 
tails 

ir?(p) ~ e -+r(p-2) f orp >3 

ir?(j>) ~ e a + T P forp<-l. 1 ) 

The exponential tails directly give f?(P) = (<t+t) _1 hxiTrip) /tw{~ p) ~ 2. Note that conversely, 
exponential tails give rise to singularities in the Fourier transform since f^dpe^Pe-M =i/(q + i)- 
i/(q-i). 

6 Discussion 

We have shown that for a dragged Brownian particle described by a Langevin equation, the work 
fluctuations satisfy the Conventional Fluctuation Theorem, and the heat fluctuations satisfy an 
Extended Fluctuation Theorem. These results may seem to only be peculiarities of this model. 
However, these results largely carry over to other more general systems too. 

1. The model of a Brownian particle in a harmonic potential dealt with here can be mapped onto 
a system of an electric circuit with a current source [10]. Furthermore, the model can be extended to 
the Rouse model for a polymer, which consists of a harmonically bound system of Brownian particles. 
It has been shown that there too the work distribution is Gaussian and satisfies the CFT [15], while 
the heat fluctuations appear to have exponential tails [16]. 

2. Another extension is to consider the non-over damped case. This has been done in the context 
of work fluctuations by Douarche et al. [17]. While the details of the work fluctuations for finite times 
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can be quite different, the general conclusion, namely, that the work fluctuations satisfy a CFT, still 
holds. It is hard to imagine that the heat would not satisfy an EFT, since this model is just a more 
detailed description of the same Brownian system. 

3. Extensions to non-harmonic potentials have been explored as well. Blickle et al. have studied 
the work fluctuations in [18], while Baiesi et al. [19] have argued for an EFT for heat under some 
general assumptions, one of which is the existence of exponential tails for the distribution of Q T . 

4. The limitation p < p* to the CFT has been found analytically by other authors as well [20,21]. 
The limit is in all these cases also due to singularities in the complex plane of the Fourier transform 
of the distribution, which correspond to exponential tails of the distribution itself. 

5. Turning to deterministic systems, Evans [22] argued that there should be a relation between 
the FT in such systems and the EFT. Indeed, Gilbert has shown numerically that the EFT holds 
in a Nose-Hoover thermostated Lorentz Gas [23] and [20] explains in a general way why the EFT 
should hold using a large deviation formalism, but still assuming that there are exponential tails and 
that the work fluctuations have faster than exponential tails. They also provide a strategy on how to 
apply the chaotic hypothesis to singular systems such as Lennard-Jones systems, by using a Poincare 
section that circumvents the singularities in the potential. 

6. A derivation of a general form of the EFT has been given for stochastic systems by Baiesi 
et al. [19] and for deterministic systems by Bonetto et al. [20], both based (among others) on the 
presence of exponential tails. They found that universal features in the EFT of the lmv^oo fr(p) 
are the slope one for p < 1 and the presence of a plateau for large p. On the other hand, the shape 
and extent of the intermediate-p region is not universal, nor is the height of the plateau, unless the 
work distribution is symmetric around its average, in which case the plateau value is always 2. 

7. Evidence for the EFT has been seen in experiments as well, notably by Gamier and Ciliberto 
in an electric circuit [7]. 

8. It has been suggested that the Gallavotti-Cohen Fluctuation Theorem (Eq. (1)), meaning 
the stationary state Fluctuation Theorem for entropy production, or heat, would not hold near 
equilibrium, but the FT for the "dissipation function" of Evans and Searles (which often coincides 
with work) would [24]. However, it seems that in the appropriate limit r — > oo, with a restriction 
p < p*, the Gallavotti-Cohen Fluctuation Theorem does hold. Nonetheless, there is a problem in 
observing the EFT for heat near equilibrium that the CFT for work does not have. As we can see in 
the exactly solved model of the Brownian particle, the FT for work decays to the correct FT form on 
a time scale of 0(r r ) [cf. Eq. ([37]) ]. One the other hand, the heat FT for p < 1 decays to the correct 
FT on a time scale of 0(a^_ 1 ) [cf. Eq. f|38j) ] . The latter time scale diverges as <r+ gets smaller, i.e., 
as one gets closer to equilibrium, whereas the time scale for the work FT (i.e., r r ) does not. This 
explains why it is very hard to see the FT for heat in simulations [25]. 
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